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We study diffusion-limited coalescence, A + A # A, in one dimension, in the presence of a 
diffusing trap. The system may be regarded as a generalization of von Smoluchowski's model for 
reaction rates, in that: (a) it includes reactions between the particles surrounding the trap, and 
(b) the trap is mobile — both considerations which render the model more physically relevant. As 
seen from the trap's frame of reference, the motion of the particles is highly correlated, because of 
the motion of the trap. An exact description of the long-time asymptotic limit is found using the 
IPDF method, and exploiting a "shielding" property of reversible coalescence that was discovered 
recently. In the case where the trap also acts as a source — giving birth to particles — the shielding 
property breaks down, but we find an "equivalence principle": Trapping and diffusion of the trap 
may be compensated by an appropriate rate of birth, such that the steady state of the system is 
identical with the equilibrium state in the absence of a trap. 



I. INTRODUCTION 

Non-equilibrium kinetics of diffusion-limited reactions has been the subject of much recent interest While 
equilibrium systems can be completely analyzed by means of standard thermodynamics methods, and reaction-limited 
processes are well described by classical rate equations , there exist no such general approaches to the problem of 
non-equilibrium, diffusion- limited reactions. 

A fundamental model for the reaction rate in diffusion-limited processes had been presented by von Smoluchowski |j| . 
In this model, an ideal spherical trap is surrounded by a swarm of Brownian particles. The rate of absorption of 
particles into the spherical trap models the reaction rate. The Smoluchowski model is limited in two important 
respects: (a) the particles react with the spherical trap but not with each other, and (b) the trap itself does not 
diffuse, but remains static at the origin. Both limitations are unphysical: In real reaction processes all particles 
interact (and may react) with each other, and all particles are mobile. 

Several attempts have been made to remove the restriction of an immobile trap |lO|-|l^] . The problem is complicated 
by the apparent correlations in the motion of the surrounding particles: With a mobile trap, the motion of the 
surrounding particles is highly correlated, since a step of the trap to the left, say, in the lab frame of reference, 
results in an apparent step to the right of all the surrounding particles in unison, in the trap's frame of reference. 
Results are restricted to empirical formulas inspired by numerical simulations Jnj, or to a number of special cases 
(immobile particles jll]]; short times jjj]]). As regards reactions between the surrounding particles, these could hardly 
be considered, other than numerically, because few models of diffusion-limited kinetics yield themselves to exact 
analysis. In fact, diffusion-limited coalescence (A + A ^ A) and annihilation (A + A — > 0) in one dimension alone 
account for most of the known exact results ]l5)-|23[]. 

Recently, we have studied reversible coalescence, A + A # A, on the line and in the presence of a static trap 
An exact analysis is possible with the method of interparticle distribution functions (IPDF) (25) . We have found a 
remarkable property of "shielding" : The particle nearest to the trap effectively shields the remaining particles from 
the trap. The steady state of the system is uniquely characterized by the distance of the nearest particle to the 
k> , trap — all other particles remain distributed exactly as in the equilibrium state of the system in the absence of a trap. 
; I ■ This shielding property persists even in the presence of a bias field (convection, or drift). 

In this paper we consider reversible coalescence with a mobile trap which diffuses with a diffusion constant Z?t, 
not necessarily equal to the diffusion constant of the surrounding particles, D. The problem may still be formulated 
with the IPDF method, in spite of the correlations induced by the motion of the trap. The shielding property of 
the nearest particle to the trap also remains in effect, and it enables us to find a complete exact description of the 
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distribution of particles in the long-time asymptotic limit. We find that relative to the trap the surrounding particles 
remain at equilibrium, but the gap to the nearest particle is proportional to D + Dt- 

An interesting generalization is to the case when the backward process A — > A + A is not limited to the particles 
alone, but the trap too may act as a source. Again, the rate of generation of particles from the trap, vt, need not 
coincide with v, the production rate from the reverse birth reaction. We find that if vt > the shielding property 
breaks down, and we are unable to derive a complete exact description of the system. However, if the birth rate is 
vt = v(l + Dt/D), the effect of the trap is nullified: The particles remain distributed as in equilibrium, as if the trap 
were not present. 

The remainder of this paper is organized as follows. In section O we review the coalescence process and the IPDF 



method used for its analysis. The model with a mobile trap is considered in section 111. In section IV, we generalize 
to the case where the trap may also act as a source. We conclude with a summary and discussion, in section [v|. 

II. COALESCENCE AND THE IPDF METHOD 

Our model |l^,[l4|,^| is defined on a one-dimensional lattice of lattice spacing a. Each site is in one of two states: 
occupied by a particle A (•), or empty (o). Particles hop randomly to the nearest neighbor site to their right or left, 
at rate D/a 2 . Thus, in the continuum limit of a — > the particles undergo diffusion with a diffusion constant D. 
A particle may give birth to an additional particle, into a nearest neighbor site, at rate v/a (on either side of the 
particle)^]. If hopping or birth occurs into a site which is already occupied, the target site remains occupied. The last 
rule means that coalescence, A + A — > A, takes place immediately upon encounter of any two particles. Thus, the 
system models the diffusion-limited reaction process 

A + A^A. (1) 

The dynamical rules of the model are illustrated in Fig. 1. 

time 
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Figure 1 Reaction rules: (a) diffusion; (b) birth; and coalescence, (c) following diffusion, and (d) following a birth event. The 
broken lines in (a) and (b) indicate alternative target sites. 

An exact treatment of the problem is possible through the method of Empty Intervals, known also as the method 
of Inter-Particle Distribution Functions (IPDF) [^5). The key concept is E n>m (t) — the probability that sites n,n + 
1, • • ■ , m are empty at time t. The probability that site n is occupied is 

Prob(site n is occupied) = Prob( • ) = 1 — E ny7l . (2) 

The event that sites n through m are empty (prob. E n>m ) consists of two cases: site m + 1 is also empty (prob. 
E n ,m+i)i or it is occupied. Thus, the probability that sites n through m are empty but site m + 1 is occupied, is 

Prob(S •••"•) = E n , m - E n>m+ i , (3) 



x Our notation here differs from previous work: we take the birth rate to be v/a rather than v/2a, to achieve a more aesthetic 
form of the final result. 
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and likewise, 

Prob(« o • • • o) = E n ,m - £ n _l, m . (4) 

With this in mind, one can write down a rate equation for the evolution of the empty interval probabilities: 

dt 
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For example, the first term on the r.h.s. of Eq. (^) accounts for the increase in E n m when the particle at the right 

edge of o • • • o • hops to the right and the sites n, . . . , m become empty; the second term denotes the decrease in 
E ni m when a particle at to + 1 hops to the left into the empty interval n, . . . , m, and so on. 

Eq. (||) is valid for to > n. The special case of m = n corresponds to E n _ n — the probability that site n is empty. It 
is described by the equation 

dE n , n _D 
_____ (l-^ tn ) 

T^(E n n E n n ^x) 
D , 

+-j i -_?„,„) 

n(En,n E n — X,n) 

— [(E n n — E n n+ i) + (E n n — E n -x „) . (6) 
a 

Comparison with Eq. (|J) yields the boundary condition: E n>n ^i = 1. The fact that the {-E^m} represent probabilities 
implies the additional condition that E n ^ m > 0. Finally, if the system is not empty then limn^-oo E n ^ m = 0. 

In many applications, it is simpler to pass to the continuum limit. We write x — na and y — ma, and replace 
E n .m(t) with E(x,y,t). Letting a — > 0, Eq. ([|) becomes 

d d 2 d 2 d d 



with the boundary conditions, 



E(x,x,t) = l, (8) 
E(x,y,t)>0, (9) 
lim E(x,y,t) = . (10) 



The concentration becomes 



d_ 

dy' 



p(x,t) = -—E(x,y,t)\ y=x , (11) 



and one can also show that the conditional joint probability for having particles at x and y but none in between, is 

P*{x,y,t) = -^-E{x,v,t). (12) 
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From P2 one obtains the "forward" (and also "backward") IPDF-— the probability that given a particle at x (y) the 
next nearest particle to its right (left) is at y (x): 

Pf(x,y,t) = p(x,ty 1 P 2 (x,y,t) ; p h (x,y,t) = p(y,t)~ 1 P 2 (x,y,t) . (13) 

The IPDF method can also handle multiple-point correlation functions fl4}| . Let E n (x\, yi, x<z, IJ2, ■ ■ ■ , %n> Vm t) be 
the joint probability that the intervals [xi, y.j] (i — 1,2, ... , n) are empty at time t. The intervals are non-overlapping, 
and ordered: x\ < y\ < ■ ■ ■ < x n < y n . Then, the n-point correlation function (the probability of finding particles at 
x\,X2, . . . ,x n at time t) is given by 

Qn 

p n (xi, ...,x n ,t) = (-1)" — — — E n (xi,yi, . . . ,x n ,y n ,t)\ yi=Xu ...,y n=Xn . (14) 

For reversible coalescence, the E n satisfy the partial differential equation: 

d d 2 d 2 d 2 d 2 

-E n (x 1 ,y 1 ,...,x n ,y n ,t)=D(^ + ^ + --- + — + —)E n 

V ^dxi <3yi + + ^dx n dy n En ' 15 

with the boundary conditions 

lim E n (xi,yi, . . . ,x n ,y n ,t) = E n -i(xi,yi, . .., ^, ii,... ,x n ,y„,t) , (16) 

XiYVi or yilxi 



and 



lim E n (xi,yi, . . .,x n ,y n ,t) = B n _ 1 (a; 1 ,j/ 1) . . .,x n ,y n ,t) . (17) 

yi1x i+1 or x i+ ilyi 



For convenience, we use the notation that crossed out arguments (e.g. ^fj) have been removed. The E n are tied 
together in an hierarchical fashion through the boundary conditions ( |l6| ) and (p~T[) : one must know £? n _i in order to 
compute E n . 

As a trivial example, consider the homogeneous steady state of reversible coalescence. This is in fact an equilibrium 
state, which satisfies detailed balance. The particles are simply distributed completely randomly — a state which 
maximizes their entropy. One obtains 

E niCq = exp{-7[(yi — x\) H h (y n - x n )]} , (18) 

and 

Pn,eq(xi,X2, ■ ■ ■ ,X n ) = 7™ , (19) 

where 7 = v/D is the particle concentration at equilibrium. 



III. COALESCENCE WITH A MOBILE TRAP 



We now consider the coalescence model but with a trap which diffuses with a diffusion constant In the discrete 
representation, the trap hops to its right or left at rate Dt/o, 2 . A particle that hops into the trap is irreversibly 
captured by it. Similarly, when the trap hops onto an occupied site it captures the particle in that site. It is 
convenient to analyze the system in the trap's frame of reference. In this view, the trap remains static at a site which 
we choose to be the origin; n — 0. When the trap does not move (in the lab frame of reference) the changes to E n . m 
are described by Eq. (|1). A motion of the trap is perceived as a coherent opposite motion of the particles in the trap's 
reference frame. Thus, the changes to E n ^ m due to the motion of the trap are: 

d D T 

"TTT (-^n,m)trap — o~ [(-^n+l.m+1 ■^n.m+l) ~l~ (-^n— l.m — 1 -^n—l.rn) 

at a A 

— (E n ,m — -EVi-l.m) ~ (E n , m — -EVi,m+l)] ■ (20) 

For example, the first term on the r.h.s. denotes the possibility that site n + 1 is occupied while the subsequent sites 
n + 2,n + 3,...,7n + l are empty, and the trap hops to the right: In the trap's frame of reference the particle at n + 1 
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seems to hop to the left, thereby clearing the [n, m] -interval. Notice that it is important to make sure that site m + 1 
is empty, since otherwise site m would become occupied as the trap moves to the right. 

Putting together all the terms in ^) and (^0|), and passing to the continuum limit, we obtain 



d d 2 d 2 



d 2 

2D T ——E 
ax ay 



V ^dx dy ^ 



(21) 



which is now valid in the infinite wedge < x < y. The term with the mixed derivative is special: it arises because 
of the correlated motion of the particles in the reference frame of the trap. 

The trap at n = could be realized by holding that site empty, at all times. Thus, it follows that E\^ m = i?o.m, 
which results in the boundary condition (in the continuum limit): 



d_ 

dx 



E(x,y,t)\ x=0 = . 



In addition, the boundary condition dlG) becomes 



lim E(x, y, t) = , 

y— >oo 



(22) 



(23) 



while (|J) and (||) still apply, without change. 

We now search for a solution of Eq. ( pl| ) , with the boundary conditions (||) , (ph , ( |22| ) , and ( p3| ) , in the long-time 
asymptotic limit, dE/dt = 0. It is simple to find eigenfunctions which obey Eq. (H), and other eigenfunctions which 
obey Eq. ( [22] ) , but we were unable to devise a systematic method for finding the linear combinations that would satisfy 
both conditions simultaneously. Instead, we offer a solution based on the newly discovered property of "shielding" in 
the coalescence model [|i 26 



In the steady state of the coalescence model with a static trap, it is found that the particles are distributed randomly, 
exactly as in the equilibrium state of the homogeneous, infinite system (end of section |l]). The system is then fully 
characterized by p{z) — the density distribution function of the distance between the trap and the nearest particle to 
the trap, z. The nearest particle effectively shields the remaining particles from the trap (Fig. 2). As we show below, 
the same shielding effect takes place even when the trap is mobile. 




Trap 



nearest 
particle 



Figure 2. Schematic illustration of the shielding effect: The particles in the shaded area are distributed randomly and 
independently from each other, as in equilibrium. The gap z between the particles and the trap follows the probability density 
distribution po(z). 

Assuming that shielding holds, let E(x,y\z) be the conditional probability that the interval [x,y] is empty, given 
that the nearest particle to the trap is at z, then: 



E(x,y\z) = { 0, 



x < y < z 
x < z < y 



e -7(»-^) ) z < x < y 

Hence, 

/■OC /"OO 

E{x,y)= / E(x, y\z)p(z) dz = / p{z) dz + e^ {y ~ x) \ p(z) dz 

J0 Jy 

= l-F(y)+e~' y< - y '^F(x) , 
where in the last equation we introduced the definition 

F(z) = [ p(z')dz' 
Jo 



(24) 



(25) 



(26) 
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If the hypothesis of shielding is correct, the particle concentration can be obtained from Eqs. ( |11[ ) and \&l 

p{x)=p{x)+ 1 F{x) . (27) 
Substituting E(x,y) from Eq. ( ^5| ) into Eq. (|2~l|), in the stationary limit, the variables separate: 

i(D + Dt)^F(x) + v^F(x)]e^ = [(D + D T )^F(y) + v^-F{y)]e™ , (28) 
and so, one is lead to the conclusion that 

{D + DT) ^ F(z) + v Tz F{z) ^ c ' e ' 1Z ' (29) 



where C is a constant. 

The general solution of Eq. (p9|) is 



F{z) = A + Be- 1 z + Ce-~> z , (30) 



where A, B, and C — C'D 2 /v 2 Dt are constants, to be determined from boundary conditions. From the definition 
of F, we have; F(0) = 0, lim^oo F(z) — 1, and F(z) > 0. The boundary condition due to the presence of the trap 
(Eq. ^) translates into dF/dz\ z= o — 0. Thus, we find 

F(z) = 1 + §-e~^ - ^^e~^ , (31) 

where 7 ' = v/(D + Z>r). It then follows that 

p(z)^(e^-e^). (32) 

From p{z) one immediately obtains the average distance between the trap and the nearest particle: {2D + Dt)/v, as 
well as the particle concentration in the trap's frame of reference (Eq. [~~ 

p(x) = 7(1 - e^' x ) . (33) 

The last result is similar to the one obtained for a static trap, only that the width of the depletion zone near the trap 
is 1/7' = (D + D T )/v, instead of I/7 = D/v. 

Our original goal of finding the empty interval probability has now been achieved. Using (^H|) and (|3l|), we get 

E{x, y) = g-^-z) + D + DT e -^' y [l - e-fr-T'K^] . (34) 

This solution can be verified by direct substitution in Eq. (|2l] ) and in the boundary conditions (||), j2^), and (|2^). 
The fact that we have found a solution proves that shielding indeed takes place, even with a mobile trap. On the 
other hand, we have merely shown that E(x 1 y) is consistent with the shielding assumption. We now wish to show 
that the same is true for the whole hierarchy of E^s, and hence for all n-point correlation functions. 

If the particles are distributed as implied by shielding, then, following a reasoning similar to that which led to 
Eq. (p5|), we should have 

E n (xi,yi, ...,x n ,y n ) = l- F(y n ) 
+e-^- x ^[F(x n ) - F(» n _i)] 
+ . . . + e -7{(»»-»«)+-+{w-*0}[ir( a;j ) _ F( w _i)] 

+ . . . + e -r{(vn-xn)+-+(m-xi)} F ( Xl ) (35) 

It is easy to confirm that these functions fulfill the boundary conditions ( fl6| ) and (|l7|) . Eq. (1^) is also satisfied, 
provided that F satisfies the same equation as above, Eq. (|29|), with the same boundary conditions. That is, the 
solution found above for F, combined with Eq. (|35|), solves the problem of the E n . Indeed, using Eqs. (p^), and 
(^), we find the rt-point correlation function: 

p n (xi, ■ ■ • ,x n ) = p(xi)-i n ~ x , (36) 

exactly as we expect from a system with the shielding property. 
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IV. COALESCENCE WITH A TRAP-SOURCE 



We now wish to consider a further generalization of the trapping problem. Suppose that the trap is imperfect, in the 
sense that it may also act as a source: the trap gives birth to A particles into the site next to it, at rate v^/a. When 
£>t = D and vt = v the trap is identical to the rest of the particles (only that the system is empty to the left of the 
trap). Thus, such a trap-source might also be viewed as a special, tagged particle, characterized perhaps by a different 
diffusing constant and back reaction rate. The model constitutes a modest first step towards the understanding of 
the more realistic situation where the size of aggregates matters, and clusters diffuse and give birth at different rates, 
determined by their accumulated mass. A very recent application is to the evolution of bacterial colonies living near 
a patch of nutrients. Nelson et al., |^7],^8| analyze such experiments with a diffusion- limited coalescence model with 
a source (modeling the nutrients) similar to ours. 

The evolution equation for empty intervals in the trap-source model is identical to that of a perfect trap (Eq. |T]) . 
The birth of particles from the trap affects only the boundary condition at x = 0: it is no longer true that the trap 
may be realized by simply holding site n = empty. To derive the appropriate boundary condition we consider the 
total changes to E n>m , which are obtained by putting together Eqs. (||) and ((2^). The case of n = 1 needs to be 
considered separately, since we do not know what is Eo m . The changes to £a, m , including birth from the trap, add 
up to: 

d D 

-Q^Ei >m — — [(Ei^m-i — Ei >m ) — (Ei. m — Ei^ m+ i) + (E2.n1 — Ei tm )\ 

D T 

^ 5"[(-^2,m+l — -Ea.m+l) — (Elm ~ E\ m +l) + (-El ,m-l — Ei m )\ 

v . . vt 

— [El m — Ex m+l) Ei m ■ (37) 

a a 

Comparison of this equation with that for general n, when n = 1, yields the discrete boundary condition: 

( 9 H ){E\ t m — E m ) H o~(-El,m-l — Eom-l) = —Eim ■ (38) 

Notice that when — this reduces to Eo jm — E\. mi as for a perfect trap. Passing to the continuum limit, the 
boundary condition for the trap-source becomes 

(D + D T )^-E(x, y, i)| a=0 = v T E(0, y, t) . (39) 



We now seek a solution to Eq. (21), in the stationary limit dE/dt = 0, and which satisfies the boundary conditions 
(|), (|), (||), and ©. Assuming that shielding holds, we follow the same steps as in section flfl and we arrive at 



exactly the same result; F(z) = A + Be 7 2 +Ce 7Z , only that now the boundary condition dF/dz\ z=0 — is replaced 
by: 

(D + D T )e^ z ^F{z)\ z=0 = v T [l - F(z)} , (40) 
az 



from Eq. (39). (Again, notice that when — one recovers the condition dF/dz\ z =o — 0.) 

From the boundary condition F(oo) = 1, we obtain A = 1. Furthermore, from F(0) — we get B = —(1 + C). 
Finally, from the boundary condition due to the trap-source, Eq. (|io|), we get 

[— 7 'e (7 - V)z - 1](1 +C) = (— 7' - l)C . (41) 
v ' v 

This condition cannot be satisfied for all z generically, and so one must conclude that shielding does not take place 
in the system with a trap-source. On the other hand, for the special case that C = —1 and ^'v^/v = 7, Eq. ( |4l| ) 
is satisfied. In this case F(z) = 1 — e~ lz , which leads to E — e~~ , ( y ~ x ^ and p(x) = 7. That is, the particles are 
distributed exactly as in equilibrium, as if there were no trap! Thus, there exists a whole class of states which are 
equivalent to the equilibrium state of the infinite, homogeneous system, without a trap. The equivalent states are 
characterized by the relation 

Dt 

v T = v(l + -^) . (42) 
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For these states the effect of the trap is nullified. A larger diffusivity of the trap, Dt, (i.e., a larger trapping efficiency) 
is exactly compensated by an increasing rate of birth vt from the trap-source. 

An interesting case is when Dt = D and vt = v. Then the trap is identical to the surrounding particles. Notice, 
however, that this is not a "special" state (Eq. [l2] is not satisfied), and hence we conclude that at the equilibrium 
state the system seems inhomogeneous from the point of view of a tagged particle. To shed some light on these 
baffling results, we first point out that since the equilibrium state is homogeneous it is perceived without change by 
any moving observer which does not interact with the particles, including random walkers. [Indeed, the equilibrium 
state E = e -7 ^ -1 -' is a steady-state solution of Eq. (|2l|).1 Imagine then an observer diffusing through the system 
with diffusion constant Dt, and which does not interact with the particles. From the point of view of the observer 
he is static, and the average concentration of particles is constant and equal to 7. Ignoring the half infinite line to 
his left, the observer could interpret crossings of particles from right to left as "trapping" events, provided that he 
also interprets crossings from left to right as "birth" events. The apparent rate of birth (crossings from left to right) 
would be v^/a = n[(D + D^)/a 2 + v/a], where n = 7a is the average number of particles at the site occupied by 
the observer. Passing to the continuum limit, we recover the "equivalence" condition, Eq. P4). (Alternatively, one 
could use the exact discrete result: 7 = v/{D + va), to obtain the discrete analogue of the equivalence condition: 
«T=«[l + i>r/(I> + t;o)].) 

Although shielding breaks down when the trap acts also as a source, one may still look for a solution to the problem 
in more conventional ways. We were unable to find an analytic solution; however, the discrete equations can be 
integrated numerically, and the particle system may also be simulated on a computer. Computer simulations confirm 
the fact that shielding breaks down when vt > 0, and that the particle concentration beyond the nearest particle to 
the trap is no longer as in equilibrium. In Fig. 3 we show typical results for for various values of vt- As vt increases, 
the concentration of particles near the trap increases, from zero (for vt = 0), to 7 (for the appropriate "special" rate, 
Eq. |4^ ), and to concentrations larger than 7. 




lattice site 

Figure 3. Concentration profile with a trap-source: Shown are results from numerical simulations (symbols) on a lattice of 
25,000 sites, averaged over a 2 /D = 10 6 time steps; as well as results obtained from numerical integration of the exact discrete 
equations (solid curves). The cases shown are for vt/v = 2, 1, 0.5, and (top to bottom), with Dt = 0. 

V. DISCUSSION 

We have studied diffusion-limited coalescence, A + A ^ A, in the presence of a diffusing trap, and we have found 
an exact description of the long-time asymptotic limit, using the IPDF method. When the trap is perfect, the system 
displays a "shielding" property: the particle nearest to the trap effectively shields the other particles from the trap. 
That is, the particles remain distributed as in the equilibrium steady-state of the infinite homogeneous system (without 
a trap), and only the distance between the trap and the nearest particle is unusual. This distance grows linearly with 
Dt — the diffusion coefficient of the trap. 

For an imperfect trap which also acts as a source the shielding property breaks down and we were unable to find an 
analytic solution, but the exact equations can then be solved numerically. We have found an intriguing "equivalence 
principle": all systems with vt = v(l + Dt/D) are equivalent to each other. The trap then seems invisible and the 
particles remain distributed as in the homogeneous equilibrium state. 

Our system is a generalization of von Smoluchowski's model for reaction rates; the particles react with each other, 
and the trap is mobile. The reaction rate equals the rate of influx of particles into the trap. This rate is k — 
n[(D + Z?t)/o 2 + v/a], where n is the average number of particles at the site adjacent to the trap. For a perfect trap, 
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we use the result of Eq. (|3^) to find n = a 2 (dp/ 'dx) x —o — a 2 ~fj' = a 2 v 2 /D(D + -Dt), and so, in the continuum limit 
(a — > 0) we get k — v 2 /D. Curiously, the trapping rate is independent of the diffusivity of the trap, Dt- A faster trap 
visits more sites per unit time, but it also depletes its immediate neighborhood more effectively, and the two effects 
cancel each other. For the case of a trap-source, we have failed to obtain an analytic expression for n, and hence we 
could compute k only numerically. 

There remain several interesting open problems. We have considered only the steady state of our model, but the 
transient is also of interest. For perfect traps the shielding property holds at all times (provided that the initial 
condition is compatible with it) and one can exploit it to find an analytic answer. An important open problem is that 
of finding a systematic method for solving the evolution equation for empty intervals. We were fortunate to come 
across a solution which obeys shielding, but shielding does not always hold, as for example for non-ideal trap-sources. 
Indeed, most remaining open questions concern the model with a trap-source. An exact analytic solution for this 
case is still missing. We have managed to prove, however, that the solution could not be of the form of a sum of 
exponentials (finite or infinite), other than for the special states, equivalent to equilibrium. 

An interesting question is whether there exist other classes of equivalence. That is, are there any states equivalent 
to each other, but not to the equilibrium state? — We have managed to prove that such states do not exist, at the 
level of the empty interval probability. However, there remains the possibility that different systems might share 
the same concentration profile, in spite of differences in their empty interval probabilities. Whether such states exist 
remains an open challenge. 
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